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ABSTRACT 

A  new  method  for  detecting  negative  cycles  in  a  graph  is 
proposed.   This  method  is  based  upon  the  primal  -  dual 
relationships  of  a  linear  program  formulated  from  an  assign- 
ment problem  type  network.   A  computer  program  is  developed 
for  this  new  method  to  include  the  complete  solution  of  the 
assignment  problem.   Results  are  given  on  program  efficiency 
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I.   OBJECTIVE 

A  new  method  has  been  proposed  by  H.  Greenberg  for  the 
detection  of  negative  cycles  in  an  assignment  type  problem, 
The  objective  of  this  thesis  was  to  develop  a  computer  pro- 
gram for  this  new  method  and  to  measure  the  efficiency 
thereof  relative  to  execution  time. 


II.   DESCRIPTION  OF  THE  METHOD 

A.   INTRODUCTION 

Suppose  there  exists  a  network  with  vertices  V  =  {l,...,n}, 
directed  edges  E  =  {(i,j)e  V  X  V},  a  non-negative  integral 
valued  function  k  giving  the  maximum  allowable  flow  k. .  over 
every  edge,  and  a  non-negative  cost  function  a  giving  the  cost 
a. .  associated  with  a  unit  of  flow  over  any  edge.   A  flow  X 
of  value  v  is  an  integral  valued  function  on  E  satisfying: 

(1)  0  <  x±.    <  kiy    (i,j)e  E. 

n 

(2)  .Z,  (x. .  -  x..)  =  V,  i  =  1, 

1=1  v  1J    J  3/    »      ' 

=  0,  i  =  2 ,  ...,  n-1, 

=  -v,  i  =  n. 

Vertices  1  and  n  are,  respectively,  called  the  flow  source 

and  sink. 

The  minimal  cost  flow  problem  is  to  find,  among  all  flows 

X  of  value  v,  one  which  minimizes: 

(3)  Q(X)  =     Z    x,,a 

(i,j)e  E   ^  ^ 

The  assignment  problem  is  to  fill  n  jobs  by  as  many  men 
at  least  total  cost  where  a. .  represents  the  cost  of  using 
man  i  (i  =  1,  . . . ,  n)  in  job  j  (j  =  n  +  1,  ...,  2n) .   Because 
of  this  bipartite  form  the  problem  can  be  thought  of  as  a 
special  minimal  cost  flow  problem. 

Now  suppose  a  flow  has  been  found  satisfying  (1)  and  (2) , 
and  for  the  assignment  problem  this  corresponds  to  finding  a 


flow  of  value  v  =  n  and  x. .  values  by,  say,  the  "Least-Cost 
Rule"  given  by  Dantzig  in  [1]  (p.  307).   We  then  construct  an 
associated  network  G(X),  given  by  Klein  [2],  which  has  the 
same  vertices  as  the  original  network  and  directed  edges  as 
follows : 

(4)  E(X)  :  (i,j)  if  xi.   <  k±.    and  x^  =  0, 

(j,i)  if  x±.   >  0, 
with  revised  capacities : 

(5)  k'  :  kr.  =  k.  .-x.  .  ,  if  x.  .  <  k.  .  and  x...  =0, 

13     iJ   iJ      iJ     13      3i 

kT •  =  x. . ,  if  x. .  >  0 , 
13     13      13 

and  with  revised  edge  costs  : 

(6)  a'  :  aT.  =  a..,  if  x. .  <  k..  and  x..  =  0, 
^  J  13     13  '     13     13      3i 

al ■    =  -a. . ,  if  x. .  >  0 . 
31     13      13 

These  revised  edge  costs  are  simply  those  associated  with 
increasing  or  cancelling  the  flow  by  one  unit  on  these  edges: 
the  revised  capacities  indicate  the  extent  to  which  this  can 
be  done. 

Now,  use  a  result  proved  in  Busacker  and  Saaty  ([3],  pp. 
256-257) . 

Theorem.   X  is  a  minimal  cost  flow  if  and  only  if  there 
is  no  directed  cycle  C  in  G(X)  such  that  the  sum  of  the  costs 
around  C's  edges  are  negative.   (A  directed  cycle  is  a 
sequence  of  distinct  directed  edges  of  the  form 

^o'1^  (il»i2^*"^ip,iq^  ^q^o^  involving  distinct 
vertices . ) 


An  immediate  consequence  of  this  theorem  is  that  a  test 
of  the  optimality  of  the  flow  X  is  at  hand  if  G(X)  can  be 
checked  for  the  existence  of  a  negative  cost  directed  cycle. 
Further,  if  such  a  cycle  is  found,  an  improved  flow  is  ob- 
tained by  sending  a  positive  unit  flow  around  this  cycle. 
Such  a  flow  alteration  obviously  leads  to  a  lower  total  cost 
and  also  leaves  the  flow  value  v  unchanged. 

B.   THE  NEW  METHOD 

Several  methods  exist  for  locating  negative  cost  directed 
cycles  (see  Klein  [2]).   The  new  method  adds  to  these.   The 
idea  behind  this  method  is  to  formulate  a  linear  program 
from  the  network  G(X)  using  the  flow  balance  equations  as 
constraints.   A  mathematical  statement  of  this  set  of 
equations  for  the  assignment  problem  is: 

(7)  minimize:   Z  =    E   x..a'.,  where  E(X)  is  deter- 

(i,j)eE  1    1-'   mined  by  the  assignment, 


subject  to: 


E    x.  .  -    E   *ij=  0,  i  =  l,---  ,n, 
j,(i,j)eE  XJ    j,(j,i)eE 


Flow 

Balance  E    x.  .  -    E  x:g  =  0 ,  j=n+l ,  •  •  •  ,  2n, 

i,(j,i)£E    1J         i,(i,j)cE 

n  2n 

Rec*Uires  J  E  E    x..    >    1, 

positive         1  .    .       ...    ij    - 


flow 


i=l   j=n+l 


x.  .  >  0. 
13  - 


Obviously,  a  feasible  solution  exists  to  this  program. 
But,  more  importantly,  if  a  negative  circuit  exists  then  the 


solution  will  be  unbounded  since  the  sum  of  the  a',  around 
the  circuit  will  be  negative  and  the  flow  around  this  circuit 
will  be  increased  with  the  solution  approaching  minus  infinity. 

The  consequence  of  this  is  that  we  have  at  hand  another 
method  for  detecting  negative  circuits  by  capitalizing  on  the 
theorem  (see  Hadley  [4])  that  an  infeasible  dual  implies  an 
unbounded  primal  when  the  primal  has  at  least  one  feasible 
solution.   Thus  the  dual  of  the  above  program  (7)  is  solved 
through  Phase  I,  only,  of  the  Simplex  Method  to  test  for 
feasibility.   If  feasible,  then  the  assignment  is  optimal. 
If  the  dual  is  infeasible,  then  the  assignment  is  not  optimal 
and  a  change  of  assignment  is  made.   This  process  is  repeated 
until  a  feasible  dual  and,  hence,  an  optimal  assignment  is 
obtained. 

The  formulation  of  the  primal  program  is  now  illustrated 
with  a  two  by  two  assignment  problem  example.   The  cost  matrix 
for  this  example  is  (8) : 

(8) 

3       4 
1 


a13 

a14 

a23 

a24 

where  i  -  1,2  denote  the  men  and  j  =  3,4  denote  the  jobs, 
hereafter  referred  to  as  assignee  and  assigned  nodes, 
respectively.   First,  an  initial  assignment  must  be  obtained, 
say  by  the  "Least-Cost  Rule".   Suppose  this  assignment  is 


(j-i),  3-1,  4-2.   Then  the  graph  G(X)  representing  this 
problem  can  be  constructed  (9) . 
(9) 

~a13X13 


a23X23 


a14X14 


a24x24 


The  flow  balance  equations  are  then  (10) : 
(10) 


Node  1 
Node  2 
Node  3 
Node  4 


X14   X13> 


'23 


24 


x13   X23' 
X24  =  x14' 


Z    "a13X13  +  a14X14  +  a23X23  "  a24X24' 


From  which  is  written  the  linear  program  (11) : 
(11)  minimize: 

subject  to: 

=  0, 
=  0 


"X13   +   X14 


23 


13 


-x 


23 


14 


24 


24 


-  o, 
=  o, 


Flow 
Balance 


X13   +   X14  +x23   +   X24  1  l* 

x.  .    >  0. 
ij  - 


Positive 
Flow 


The  dual  to  this  program  (11)  is  then  solved  through  Phase  I 
to  test  for  optimality  of  the  assignment. 
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III.   PROGRAMMING  THE  ALGORITHM 

A.   THE  STANDARD  FORMAT 

The  dual  program  to  (11)  is  (12) : 

(12)  maximize  v., 
subject  to: 

+  v  < 

~a13' 
£  a14> 
-  a23> 
±   "a24> 

which,  after  conversion  to  the  standard  form  for  solution 
is  (13) : 

(13)  maximize  v^, 


vl 

+v3 

vl 

"v4 

v2 

'v3 

"v2 

+v4 

subject 

to 

: 

vi 

'V3 

"sl 

= 

a13' 

vi 

"v4 

s2 

= 

a14' 

v2 

'V3 

'  S3 

= 

a23' 

v2 

"v4 

"s4 

- 

a24J 

with  slack  variables  s..  through  s..   Examining  the  array  of 
coefficients,  a  standard  format  for  any  size  problem  can  be 
detected.   This  format  can  be  illustrated  with  a  three  by 
three  problem  annotated  with  general  dimensions  (14)  ,  where 

n  denotes  the  number  of  assignee  or  assigned  nodes  (in  this 

2 
case  n  =  3),  v.  (i=l , • • • ,2n)  the  dual  variables,  s.(j=l,,,,,n  ) 

the  slack  variables,  -s.  the  artificial  variables  added  to 

3 


11 


(14) 


Basic 
Variables 


n 


V1V2V3 


■n 


V4V5V6 


<- 


n 


-=>■ 


S1S2S3S4S5S6S7S8S9 


Constants 


A  f 


n 


I 


-s 


-1 


14 


L15 


16 


f 


n2 


1 


-1 


-1 


'24 


25 


26 


V 


♦ 


n 


I 


-s, 


-1 


-1 


34 


35 


d36 

■Ca15  +  a24  +  a36} 


•w 


1-1-1 


111 


0  10100001 


obtain  a  basis,  and  W  the  infeasibility  form.   Examining  the 
format  in  detail,  it  is  noted  that  the  first  n  columns  derive 
from  the  flow  balance  constraints  on  the  assignee  nodes 
(numbered  from  1  to  n) ;  the  second  n  columns  derive  from  the 
constraints  on  the  assigned  nodes  (numbered  from  n+1  to  2n) ; 
and  the  n   slack  columns  reflect  the  assignment.   The  slack 
columns  reflect  the  assignment  in  that  the  rows  in  which  the 
-l's  appear  indicate  the  assignment  that  has  been  made.   If 
the  rows  in  each  group  of  n  rows  are  numbered  identically  from 
n+1  to  2n  then  the  row  number  in  which  the  -1  appears  in  the 
first  group  of  n  rows  is  the  number  of  the  node  assigned  to 
the  first  assignee  node.   Similarly,  the  remainder  of  the 
assignment  can  be  determined.   In  this  example  (14),  5  is 
assigned  to  1 ,  4  is  assigned  to  2,  and  6  is  assigned  to  3. 
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Of  course,  only  n  of  these  -l's  may  appear.  Moreover,  the  row 
number  of  a  -1  may  not  be  duplicated  between  groups  since  this 
would  indicate  a  duplicate  assignment. 

It  is  seen  then  that  the  only  format  change  between  dif- 
ferent problems  is  in  the  slack  columns  reflecting  the 
assignment.   This  feature  was  exploited  in  programming  the 

simplex  iterations  by  reducing  the  computer  storage  require- 

2 
ment  for  this  matrix  to  a  single  vector  of  dimensions  n  by 

1.   This  vector  records  the  assignment  and  its  contents  are 

scanned  when  format  information  is  needed.   It  should  also  be 

noted  that  with  this  format  the  matrix  of  coefficients  .is 

unimodular.   Hence,  for  every  iteration  the  coefficients 

will  be  +  1  or  0 . 

B.   ITERATIVE  PROCEDURE  USED 

The  revised  simplex  method  was  used  to  perform  the  Phase 
I  iterations  because  of  the  computational  advantages.   Exam- 
ining the  standard  format  (14)  ,  we  see  that  at  least  n 
iterations  are  required  because  of  the  -l's  in  the  first  n 
rows  of  the  infeasibility  form  equation.   To  reduce  compu- 
tational time  the  revised  simplex  method  was  modified  to  carry 
out  these  n  iterations.   The  simplex  multipliers  were  not  used 
to  find  the  candidate  for  the  new  basic  variable.   Instead,  the 
first  n  variables  were  introduced  into  the  basis  sequentially. 
Subsequently,  the  standard  revised  simplex  method  was  used. 

At  this  point  it  would  be  useful  to  illustrate  the  contents 
of  the  matrix  used  to  perform  the  iterations.   Again,  this  will 
be  illustrated  with  the  previous  three  by  three  problem,  anno- 
tated with  general  dimensions  (15)  .   This  matrix  is  denoted  B 

13 


(15) 


<^ 


n 


Sl"S2S3"S4S5S6S7S8"S9 


Constants 


ACOL 


VCOL 


A 


n 


M 


14 


15 


16 


24 


25 


26 


34 


35 


36 


7 

-8 

9 

-10 

11 

12 

13 

14 

-15 


0   -10-1000   0-1 


(a15+a24+a34) 


in  the  computer  program.   In  addition  to  containing  the  in- 
verse of  the  basis  and  the  simplex  multipliers,  where  -s. 
denotes  the  artificial  variable  added  for  the  corresponding 
negative  slack  variable,  the  matrix  contains  an  assignment 
column  (ACOL)  and  a  variable  identification  column  (VCOL). 
The  assignment  column,  with  binary  entries  of  0  or  1  only, 
where  l's  record  the  assignment  (in  this  case  the  initial 
assignment  is  1-5,  2-4,  3-6),  is  used  as  a  memory  device  in 
performing  the  iterations  and  in  printing  the  solution.   The 
variable  column  records  the  variables  in  the  basis  and  is  up- 
dated after  each  iteration.   The  minus  indicates  an  artificial 
variable.   If  a  problem  is  infeasible,  this  column  is  scanned 
to  determine  which  artificial  variables  remain  in  the  basis. 
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This  information  is  necessary  in  changing  the  assignment, 
which  shall  be  described  later. 

While  the  revised  simplex  method  requires  less  computer 
storage  than  the  standard  simplex  method,  the  storage  require- 
ments for  this  method  are,  nonetheless,  relatively  large.   For 
instance  a  15  by  15  problem  requires  approximately  210,000 
bytes  of  storage.   (228  X  227  =  51,756  storage  locations  for 
the  matrix  B  +  675  storage  locations  for  a  computing  register 
and  input/output  requirements  =  52,431  total  X  4  bytes  per 
location  for  single  precision  =  209,724  bytes.)   But,  other 
methods  for  detecting  negative  circuits  require  as  much 
storage . 

C.   INITIAL  ASSIGNMENT 

Up  to  this  point  an  initial  assignment  has  been  assumed, 
but  the  method  for  finding  that  assignment  has  not  been 
described.   The  method  selected  was  the  "Least-Cost  Rule" 
with  steps  as  follows : 

1)  The  cost  matrix  is  scanned  to  determine  the  lowest 
cost,  i.e.,  find  Min  a... 

2)  The  column  j  and  the  row  i  in  which  this  cost 
appears  constitutes  the  first  assignment.   That  column  and 
row  is  then  ineligible  for  further  assignment. 

3)  Steps  1)  and  2)  are  repeated  for  costs  remaining 
in  eligible  rows  and  columns  until  all  assignments  have  been 
made. 


15 


D.   CHANGE  OF  ASSIGNMENT  RULE 

Summarizing  the  procedure  described  thus  far  -  First,  an 
initial  assignment  is  determined.   Next  the  dual  problem  to 
(7)  is  solved  through  Phase  I  of  the  simplex  method.   If  the 
problem  is  feasible,  then  the  initial  assignment  is  optimal 
and  the  solution  is  printed.   If  the  problem  is  not  feasible 
then  the  initial  assignment  is  not  optimal  and  a  change  in 
assignment  is  required. 

The  key  to  devising  an  assignment  change  rule  is  contained 

in  the  artificial  vectors  remaining  in  the  basis  because  they 

2 
reflect  the  infeasibility .   In  every  column  of  the  last  n 

columns  in  which  an  entry  other  than  zero  appears  in  the 

artificial  vector,  the  sign  is  changed  of  the  +  1  in  that 

column  of  the  basic  format  (14).   Thus,  the  artificial  vectors 

remaining  in  the  basis  are  identified  and  constructed  from  the 

inverse  and  the  procedure  just  described  is  performed  for  each 

of  these  vectors  sequentially.   The  efficacy  of  this  rule  is 

evaluated  in  the  following  section.   It  should  be  noted  that 

since  this  procedure  does  not  change  the  coefficients  of  any 

variable  in  the  basis,  the  inverse  for  the  last  tableau  remains 

valid.   Therefore,  it  is  not  necessary  to  cycle  back  to  the 

initial  tableau.   Rather,  after  the  assignment  has  been  recorded 

in  the  memory  device  (ACOL) ,  the  simplex  iterations  can  continue 

from  the  last  tableau. 

A  good  illustration  of  application  of  the  assignment  rule 

is  problem  5  (see  Appendix  B) ,  a  five  by  five  problem.   After 

nine  iterations  the  program  was  found  to  be  infeasible  and  two 
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artificial  vectors  remained  in  the  basis  (-s,R  and  -s,.). 
The  vectors  shown  on  the  following  page  trace  the  change  in 
assignment  from  initial  contents  of  ACOL  to  final  contents 
of  ACOL.   For  ease  in  recognizing  the  assignments,  the  vectors 
have  been  partitioned  by  n  (n  =  5) .   Recall  that  the  entries 
of  ACOL  are  binary  (0,1)  and  that  an  entry  other  than  0  in 
the  artificial  vector  "flip-flops"  the  corresponding  ACOL 
entry. 
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IV.   RESULTS 

Twenty  problems,  ranging  in  size  from  four  by  four  to 
15  by  15,  were  tried.   The  entires  comprising  the  cost  matrix 
were  the  integers  one  through  nine,  selected  randomly  (see 
Appendix  B) .   The  method  yielded  an  optimal  solution  in  each 
case.   The  execution  times  ranged  from  0.7  seconds  for  a 
four  by  four  to  56  seconds  for  a  15  by  15  (see  Appendix  A) . 
In  three  instances  the  initial  solution  was  optimal  and, 
hence,  an  assignment  change  was  not  required.   But,  in  the 
remaining  17  problems  an  assignment  change  was  required  and 
in  every  case  the  change  rule  produced  the  optimal  assignment 
on  the  first  attempt.   With  respect  to  the  number  of  iter- 
ations required  to  detect  a  negative  circuit  the  range  was 
from  7  to  33.   And,  in  all  cases  the  number  of  iterations 
required  was  less  than  the  number  equations  involved. 
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APPENDIX   B 
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ATQTX,  NOTE  T 
A  15  X  15  PRO 
0)  NPROB 

MEASURE  EXEC 

ME 

1»50 


?30) , L,&rnL,VCOL,K,M, 
L(230) ,M230,230)  ,0 

THF  nATA/ZERH  the  INVERSE  AND 
GNM.PNT.  THE  DATA  ARE: 
M=NUMR«=P  n«=  ASSIGNED  NnnES, 
HIS  PROGRAM  IS  DIMENSIONED  TO 
RLEM. 

UTION  TIMES. 


THE  ASSIGNMENT  FOR  PRINTING  PURPOSES  AND  TO 

1,?30 

1,230 
0 

0)  N 
,N 

0)  (X(I ,J) ,1=1, N) 
0) 
TX  (X,N,N) 

HfiTA. 
0 
S  NUMBER  OF  ASSIGNMENTS  MADE. 

Q 

i  ,n 

1,N 

(K)-N)  GO  TO  4 


1  ,N 

T.O)  KO    TO  3 

.GE.MIN)  GO  T 

,J) 


0  3 


JCOL    +    N 
ICOUNT    +    1 
.LT.N)    GO    TO 
^0) 

I  ,N 

50)     I,A(T) 

INES    LOAD    THE 

DE    WHEN    EACH 

1 

NUMBER    OF    ITERATIONS. 


N+? 
N+3 


INVERT  HE  THC 
ELEMENT  OF  B  IS 


RASTS  (B) 
LOADED. 


L. 


=  1 


If" 


PO 
J  = 
I  = 
)  - 


11 


W  EXCFDT  THF 

1,N 

1  ,N 

N.EO.  I)  GO 
0 


CONSTANT  (SIMPLEX  MULTIPLIERS). 


TO  9 


=  -1 


+  1 


29 


c 
c 
c 


c 
c 
c 


c 
c 


13 

14 


LOAD  COLU 
K  =  1 
J  =  1 
00  12 
B(I  tU 
J  =  J  + 
IFC J.L 
K  =  K  + 
J  =  1 
12   CONTIN 

LOAD  ASSI 
DO  14 
J  =  (K 
DO  13 
TF( A(K 
B( J  +  I  , 
GO  TO 
CONTIN 
CONTTN 

LOAD  CONS 
DO  16 
IF( B( T 
BUtVC 
GO  to 
B(  I,  VC 
CONTJM 
SUM  = 
DO  IP 
1MPC I 
GO  TO 
SUM  = 
CONTJM 

B(L,L) 
NCXT  2 
^IRST 
MINI  Ml) 
M  =  1 
K  =  1 
KSTAP 
MIN  = 
DO  21 
IP(B( I 
MIN  = 
PROW  = 
CONttm 

UDD*TF  VC 
B(PROW 
&ERFORM  I 
DO  81 
DO  23 
IF( I  .F 
R( I ,J) 
GO  TO 
B(I,  J) 
CONTIN 
BCL, J) 
CONTIN 
IN 
M  + 
K  =  K  + 
KSTAR 
T  F  (  K  .  L 
NEXT  2 
TNFFASIBIL 
COUMTFR. 


MN  OF  CONSTANTS  EXCE°T  LAST  ROW. 


17 
18 

THE 
THE 
THE 


19 


21 


22 
23 

81 


UPDATE 
M  = 


THE 


4° 


CM  IN  = 
JCONT 
JST*R 
DO  2? 
CRAR  = 
NEXT  4  LI 
DEPENDING 


I  =  1,0  •■:'•■"    . 

=  X(KtJ) 

1  ■::  . 

F.N)  GO  to  12  •.->••■ 

1  .  ■,.  .::     '  ,_.  ..  '. 

U^  .  -\ 

GNMFNT  COLUMN  (ACOLJ . 

K  =  1,N 

*N)  -N 

I  =  1,N  .  • 

)  -  N.NE.I)  GO  TQ  13 

ACOL)  =  1 

14 

Uc 

UF  •     .   N 

TANT  COR  LAST  ROW.  • 

T  =  1,0 

tACOD.EO.l)  GO  TO  15 

OL)  =  2*N  +1 

16 

OL)  =  (2*N  +  I )  *(-l) 

uf  :  .. 

0  .•' 

1  =  1,0 

,ACOL)  .^0.1)  GO  T0  17 
18  ■ 

SUM  +  8(1  ,L) 
Uc 
=  com  *  (-1) 

3  LIMES  M&KF  TH=  «=TRST  M  S!MD|_EX  ITFRATTONS. 
0  LTNFS  FIND  THC  PIVOT  pnw  (P<?HW)  RY  FINDING 
M  CONSTANT  VALUC  IN  EACH  GROUP  OF  N. 


=  K  +  (N-l) 

QQ3Q 

i  '  =  'ktkstar 

iD.GE.MIN)  GO  TO  21 
PC  I  ,L) 

I  -.-  '"  " 

UF 
OL. 

tVCOL)  =  M 
TFPATION. 
J  =  1,L 

I  =  K, KSTAR         .  ....-■ 

0.DROW)  GO  Tn  22 

=  PCI ,J)  -  BC  DRHW,J) 
23 

=  BC1,J) 
UE 

=  BCLtJ)  +  PCPROW,J) 
UE 

DEXES. 
1 
N 

=  K  +  (N-l ) 
T.N*N)  GO  Tn  19 

4  LTNCS  FIND  THC  MINIMUM  VALUC  IN  THF 
Ty  FORM  EOUATION  CCMIN).   JCONT  IS  A  COLUMN 
JSTAQ  IS  AN  INHEX. 

1 
=    1  -■- 

=    D    +    2*N 
J    =    1,JSTAR  •--  ". 


0 

NES 


DETFDMTNE    THE    ARITHMETIC    Tn    RE    DONE 


UPON    LOCATION    TN    STANDARD    FORMAT. 


30 


IF(J.LC.N)  GO  TO 
IF(J.tP.2*N)  CO 

IF(B(  JCGN-rt&CGL) 


24 

Tfl  ^ 

.EO.'l)  GO  TO  26 
+  B(L,JCONJT) 


B(L,I+K) 


N 
B(L,T+K)*(-1 ) 


GO  TO  29 


ST  Tp  o£TFPMTN<E  TP  AN  ITERATION  !S 
MG<=  TM  ASSIGNMENT  is  REOUIRpn.  I  <=  not, 
TEN  AND  THE  PROBLEM  rOUNTEP  IS  TESTED. 

TO  31 
GO  to  43 


CBAR  =  CBAR 
GO  TP  27 

24  K  =  <J-1)*N 
no  cl  I  =  1,N 

51   CBAR  =  CBAR  + 
no  to  ?8 

25  K    =    J    -    (N+l) 
DO    53     I    =    1,0, 

53       CBAR    =    CBAR    + 
GH    tp    2« 

26  C*AR    =    CBAP    +B(L,JCONT)    *(-l) 
?7       JCOMT    =    JCONT    + 
28       TF(  CBAR.GF  XMI  M) 

CMTN    =    CBAR 

JCOL  =  J 
2°   CONTINUE 
C  THE  NEXT  21  LINCS  TF 
C  REOUIREO  OP  TF  A  CHA 
C  THE  SOLUTION  IS  WRIT 

ITFR  =  ITFR  +1 

IF(CMIM.LT.O)  HO 

IP(BCL»L) .NE.O) 

WRITE (6,60)  ITER 

K  =  1 

DO  86  I  =  1 

DO  R7  J  =  1 

TF  (BU+J-1 

A(K)  =  N  +  J 

K  =  K  +  1 
87   CONTINUE 
86   CONTINUE 

DO  88  I  a  1,N 
8«   WRITE  (6,50)  I,A(T) 
C   MEASURE  ttmp  frhm  DEtccttoN  Oc  NEGATIVE  CIRCUIT  UNTIL 
C   SOLUTION  OF  opobl^M, 

CALL  RETIME  (IPT 

EL  =  IET  *  .0000; 

WRTTP  (6,120)  ND| 

WRITE  (6,130)  EL 

NPROR  =  NDROP  - 

IF  (NPROB.GT.O) 

STOP 

NEXT  2°  LTMES  COr 

ARS  (PCOL)  AND  c 

IF< JCOL.LE.N)  GQ 

IF( JC0L.LF.2*N) 

K  =  JCOL  -  ?*N 

Ic(B(Kt ACOL) .EQ 


,N 


n 
N 
ACOL) .EO.O)  GO  to  87 


) 

26 

ROB 


no  TO  85 


C  THE 
C  APPE 
31 


INSTRUCT  THE  COLUMN  IN  WHICH  CMIN 
INOS  THC  PIVOT  ROW  (PROW). 

TO  34 
HO  TO  35 


1,0 

B(I 

1,0 
B(  I 


i,n 


DO  37  I  = 

37  PCOL (I )  = 
GO  TO  36 

32  On  33  I  = 

33  PCOL  (I  )  = 
GO  Tn  36 

34  K  =  (JCOL-1) 
DO  3"  I  = 
PCOLU  )  = 
DO  3  8  J  = 

38  PCOLU)  = 
GO  TP  3  6 

35  K  =  JCOL  -(N+l) 
DO  3°  I  =  1,0 
PCOL (I  )  =  0 

DO  3  °  J  =  1 , 0  ,  M  ■ 
30   PCOL (  I  )  =  PCOL (I 

36  MIN  =  9999 

no  M  I  =  1,0 
IF(PCOL(I ) .PO.l) 
GO  TO  41 


1)  GO  to  32 


K) 


K)*(-l) 


0 

1,N 

PCOLU 


)  +  R(T,K+J) 

)  -  R(!,K+J) 
GO  T0  42 


31 


42 


THE 
41 


C 

c 


c 
c 


72 


IF(B(I ,L).GE.MIN)  GO  TO  41 

MIN  =  B(  I  ,U 

PROW  =  1 

NEXT  17  LINES  PERFORM  THE  SIMPLEX  ITERATION. 

CONTINUF 

B(PROW,VCOL)  =  JCOL 

DO  76  J  =  1»L 

DO  75  I  =  1,D 

IF( I  .EQ.PROW)  GO  to  72 

IF(PCOL(I ) .EQ.-l )  GO  TO  73 

IFCPCOLU  )  .FO.O)  GO  TO  72 

IF(PCOL( I ) .FO.l)  GO  xO  74  . 

GO  to  75 

B(I,J) 


73 

74 

75 

76 

THE 


75 


=  B(I  ,J)  + 


75 


B(I,J) 
8(L,J) 


B( I  ,J) 

GO  TO 

B( I ,J) 

GO  TO 

B(I»J)  = 

CONTINUE 

BU,J)  = 

CONTINUF 

GO  TO  49 

N^XT  2  5  LINES  MAKES 
GETIMP  MEASURES  TIME  TO 
43   CALL  GETIME  (IFT) 

EL  =  IET  *  .000026 

WRITE  (6,110)  EL 
RESET  TIMER. 

CALL  SETIMP 

DO  48  I  =  ltD 

IF(B(I ,VCOL) .GP.O) 

DO  47  J  =  1,D 

IF(B( J,ACOL) .FO.l) 

AROW(J)  =  B(!,J) 


B(PROW,J) 
B(PROW,J) 
(CMTN*B(  PROW, J)) 


THE  CHANGE  IN  ASSIGMPNT. 
DETECT  NEGATIVE  CIRCUIT. 


GO  to  ^b 
GO  TO  44 


to  46 


GO 
GO 


TO 
TO 


83 
84 


=  0 


=  1 


AROW  (J)  STORES  CONTENTS 
ARTIFICIAL  VECTOR. 
GO  Tn  45 

APHW(J)    =    8(1,  J )*(-!) 
IF( AROW( J) .NE.O)        GO 
GO    to    47 

IF(B( J,ACOL) .FO.l) 
I«=(B(J,ACOL) .EO.O) 
GO    to    47 
B( J, ACOL) 
GO    TO    47 
B( J,ACOL) 
CONTINUc 
WRITE( 6, BO) 
WRITE     (6,90) 
CONTINUE 
WRITE<6,70) 
GO    TO    49 
FORMAT(  HO) 
FORMAT     (1018) 

(1H1,10X, 

( ////) 

(1H  ,14, T6,  «-',T7tTA) 

(lH0t5X 

•OPTIMUM 

(1H0,5X 


OF  LAST  M*N  COLUMNS  OF 


44 
45 

46 


83 

84 
47 


48 


10 
20 

30 
40 
50 
60 


8( I ,VCOL) 
(APOW(K) 

ITER 


K  =  1,0) 


70 


FHRMAT 
FORMAT 
PORMAT 
FORMAT 
114, 5X, 
FORMAT 


•INPUT  DATA') 


llTERATION-',!^ 


•  PROBLEM  FEAMPLF 

&SSTRNMENT  I ^ • ' ) 

•PROGRAM  TNFEASIBLE 


AFTER  ITERATION-' , 


AFTFP 


80  FORMAT 
90  FORMAT 
100  cOPMAT 
110    FnR^&T 


(1H0,5X 
(1H0,5X 
(]HO,*X 
( 1H0,5X 


'ARTIFICIAL 

40T3) 


VECTOR    FOR    -ST4) 


1WAS« ,E16.6,'SEC 


•COST 
«TI  MP 
») 

•TIME 


IS 

TO 


SIM 

DETECT 


NEGATIVE  CIRCUIT 


120  FORMAT  (1H0,5X, 
1 'POLLOWS1) 
1**0   FORMAT  (  1H0,5X,F16.6, 'SEC  ) 
END 


TO  SOLVE  PROBLEM-' , 14, 5X 


32 


SUBROUTINE  WRTMTX  ( K f T ROW, JCHL ) 
'DIMENSION  K(50,50) 
WRITF  (6,30)  ( I ,  I  =1 tJCOL) 

30  P0RM&T  (//  1H  ,6-X,  M/J  =  «  t50T3) 
WQITC  (6,31) 

31  FORMAT  (1H  ) 

32  c0RMAT  (1H  ,3X,T2,3X,«50I3) 
On  33  T  =  1,  TRHW 

33  WRITE  (6,32)  I,  (K(I,J),  J  =  lfJCOL) 
RETURN 


33 
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